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In the dual formulation of d=3 SU(2) LGT, the link variables are group representations and valid configurations 
are those satisfying a number of triangle inequalities. In algorithms for local updates that automatically respect 
these constraints were described. It was also pointed out there that these local updates were not ergodic. In this 
presentation, we describe two different quasi-local updating algorithms which, in conjunction with the local 
updates, appear to be ergodic. 



1. THE DUAL FORMULATION 

A brief introduction to the dual formalism for 
d = 3 SU{2) lattice gauge theory as well as 
to the techniques that have been developed for 
its numerical simulations was presented at LAT- 
TICE'99 §. As was stated there, the partition 
function of the conventional LGT can be con- 
verted, upon using the techniques of character 
expansion and group integrations, into the dual 
form 1^ 

This partition sum is defined over the dual lat- 
tice where {ja} live over the links and {jb} over 
the diagonals to the plaquettes. The conven- 
tion for the diagonals is that they connect the 
vertices of the odd sublattice. We have collec- 
tively designated {ja},{jb} by {j}. Each cube 
of the dual lattice is spanned by 5 tetrahedra 
of which one is spanned entirely by {jf,} while 
four are spanned by three {ja} and three {jb}- 
In what follows, the b-links will be shown by 
dashed lines while the a-links will be shown 
by solid lines. Each tetrahedron carries a weight 

factor which is the SU (2) 6-j symbol | ^ ^ ^ j 

Periodic b.c for the original lattice is crucial for 
this construction. 
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2. QUASI-LOCAL UPDATES 

An important ingredient in the techniques for 
numerical simulation of such classes of models 
is an updating technique that respects all the 
triangle inequalities. In |jl[| we introduced two 
distinct classes of updating algorithms: local 
and quasi-local. The local algorithms work 
very much like the usual updating techniques for 
Monte Carlo simulations: all links but one are 
held fixed and the free link is updated by either 
Metropolis or heat bath methods. In the present 
context, one can see immediately that such lo- 
cal updates are not ergodic. To see this, let us 
recall that the representations of SU{2) can be 
grouped into two classes i.e half-integral which 
can be called fermionic and integral which can 
be called bosonic. Let us call a move that takes 
fermionic links to bosonic, and vice versa, a Z2- 
flip. It is obvious that in every triangle exactly 
two links must be Z2-fiipped simulateneously if 
we are to maintain the triangle inequalities. As 
the local moves change only one link at a time, 
they can not accommadate Z2-fiips and the moves 
therefore never take configurations out of their 
Z2-classes. 

To circumvent this the quasi-local moves were 
introduced. Since in each triangle exactly two 
Z2-flips have to be carried out, these flips can 
proliferate throughout the lattice unless a way out 
is found. Therefore one can look for those clusters 
of 2'2-flips that have the smallest volume possible. 
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In this way we identified the globeven moves at 
even-sublattice points which flip the six interior a- 
hnks as shown in fig.l. All the "exterior" b-links 
are held fixed during a globeven move and this 
contains the Z2-flips to the interior. It should be 
noted that eight tetrahedra get updated during 
each globeven move. 
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3. KAGOME VARIABLES 

Though we have identified the geometry of the 
cluster of links to be updated during quasi-local 
moves, we have to address the issue of how the 
new spin values are to be chosen such that all tri- 
angle inequalities are respected. Let us first recall 
the form of triangle inequalities. Considering the 
triangle (ji , j2, ji2) of fig. 3, these inequalities take 
the form 

\jl-j2\ < jl2 < jl+j2 
\jl - jl2 1 < 32 < jl + jl2 

\j2-jl2\ < jl < j2+jl2 (2) 



Figure 1. Quasi- local update at even sites. 



While the globeven moves are capable of Z2- 
flipping the a- links, they leave the b-links un- 
changed. For this reason, these moves are also not 
ergodic. Since b-links connect only odd sites, any 
quasi-local move affecting them must take place 
at the odd sites. It then follows that the smallest 
volume of Z2-flippings involves all the 6 a-links 
and 12 b-links at the odd site while keeping all 
the other a-links and b-links fixed with 32 tetra- 
hedra getting updated. This is illustrated in fig. 2. 




Figure 2. Quasi-local updates at odd sites. 



Figure 3. A Single Tetrahedron. 

We illustrate the K-variable strategy with the 
example of a single tetrahedra. The idea behind 
these variables is to replace inequalities by equal- 
ities to some extent. They were first introduced 
by Bargmanjp| and were later developed by An- 
ishetty et al |J. Again, let us first introduce them 
for the triangle considered above. Their construc- 
tion involves a pair (ni,n2) for each link or j- 
value such that 

«i + "-2 = jl 
«i + n-3 = j2 

nz + n2 = ji2 (3) 
These relations are invertible: 

2"! = jl + j2 - jl2 
2^2 = jl2 + jl ~ j2 

2^3 = j2 + jl2 - jl (4) 
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From eqn (4) it is clear that rii are either integers 
or half-integers and also positive semi- definite. 
The moral of the story is that wc can choose any 
triplet (711,712,713) and the ji constructed from 
them according to eqn (3) will automatically sat- 
isfy the triangle inequalities of eqn (2). 

Now let us consider the tetrahedron of fig. 3 and 
ask how we can choose ( ji , j2 , is ,ji2, j23 , jis ) such 
that the four sets of triangle inequalities are all 
satisfied. At the level of ji this can be done as 
follows: freely choose 3 j's originating at a ver- 
tex eg {ji, ja)- Next choose ji2 in the range 
given in eqn (2) and choose J23 likewise. Finally 
choose ji3 such that it simultaneously satisfies the 
triangle inequalities relevant for the two triangles 
(ji, is, jis) and (ji2, J23, jia)- This procedure gets 
increasingly tedious. Not only is this procedure 
tedious, it docs not even guarantee a solution al- 
ways even though in the case of the single tetra- 
hedron a solution is always possible. 

That a solution using this strategy may not al- 
ways be possible can be seen by applying it to 
the gobeven geometry of fig.l. Here one could 
start by first selecting a Z2-flippcd value for, say, 
OE. The link OA could then be chosen to satisfy 
the inequalities for OAE, the link OB chosen to 
satisfy OBE and OAB. So far the situation is like 
that for a single tetrahedron and all would be fine. 
Next, OC will have to be chosen to satisfy OCE 
and OBC; already at this stage there is no guar- 
antee that these two sets can be simultaneously 
realised. OD will have to be chosen to satisiy 
ODE, AOD and COD; OF will have to satisfy 
FOA, FOB, FOG and FOD! The situation only 
gets worse with these additional requirements. 

Now we show how the same problem can be effi- 
ciently addressed by using the K-variables. Again 
we illustrate the method by first applying it to 
the simple case of a single tetrahedron. For this 
purpose we introduce the variables (711,712,773) 
for the set (ji,j2,Ji2), (714,775,776) for (j2,j3,j23), 
(777,773,779) for (ji,.i3,.ii3), and (h^h^h) for 
(ji2j J23, Jis)- Writing down the various equali- 
ties we get 

jl = 712 + 713 = 778 + 779 
j2 = 711+713 = 775 + 776 

is = 714 + 716 = 7I7 + Tig (5) 



and 








ii2 


= 771 + 772 


= l2 


+ h 


i2S 


= 774+775 


= h 


+ I3 


iis 


= 717 + 778 


= h 


+ h 



Prom eqn(5) it follows that (771,772,713,774) can 
be chosen to be unrestricted, one of (775,776) can 
be chosen to be independent but with restricted 
range and one of (717,718,779) to be also inde- 
pendent but with restricted range. Choosing 
(716,779) to be the independent ones, the relevant 
restricted ranges are 

< 7l6 < 771 + 773 

< 7I9 < 77im(7l2 + 713,714 + 716) (7) 

Though (^1,^2,^3) are determined by the rii, the 
fact they all have to be positive semidefinite yields 
additional restrictions which are best seen by 
solving for Zj using eqn (6): 

h = 7l3 + 774 - 779 
I2 = 7l2 + 7l6 - 7I9 

h = Til + 7I9 - 776 (8) 

Summarising, the solution to the single tetra- 
hedron problem in terms of the K-variables is: 
(711,712,713,714) are unrestricted except for pos- 
itive semi-definiteness, (716,719) are independent 
but restricted in range and (715,777,778) are de- 
pendent: 

< 716 < 771+713 

716 — 77i < 779 < min{n2, 774) + min{n3, uq) (9) 
and 

775 = 77i + 773 - 776 

717 = 7I4 + 776 — 

718 = 7l2 + 77s - ng (10) 

3.1. TWO TETRAHEDRA CELL 

Our solution to the single tetrahedron given by 
eqn(9) will form the basis for the eventual solu- 
tion of the globeven and globodd problems. The 
final solution will be arrived at in three stages: 
i) the solution of the two-tetrahedra problem, ii) 
the solution of the four-tetrahedra problem and 
iii) globeven and globodd using the results of ii). 
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We now show how the two-tetrahedra problem is 
solved. For this purpose, consider two tetrahe- 
dra glued together as shown in fig. 4. We take 
the common triangle ABC with AB = ji, BC = 
j2,CA = j\2 to be associated with the inde- 
pendent K-variables (711,712,713). For the second 
tetrahedron ABCD we have to introduce 6 more 
variables but there will be 3 more relations as- 
sociated with the links AB,BC,CA thus reducing 
the additional variables to 3. 



Figure 4. A two-tetrahedra cell. 



The spin assignments arc taken to be : 
BE = j3,AE = ji3,CE = j23 and BD = 
k3,AD = k\3,CD = k23 and the additional 
K-variable assignments to be : {n^, 715,715) 
with (j3,j2,j23), {n7,ns,nQ) with (ji,j3,ji3), 
{1714,1715, me) with (^3,^2,^23), and {mr,ms,mg) 
with {ji,k3,ki3). Wc give only the final results: 
the unrestricted variables are (711,712,113,715,7715) 
and the restricted but independent variables are 
(rie, rig, 7716, TOg) and the remaining dependent 
variables are {n4,7iT,7is,m4,mT,ms). The equa- 
tion corresponding to eqn (9) are: 

< ne < ni -|- ns 
fiQ — Til < rig < min{n2, 175) + min{n3, rig) 

< rriQ < 711+713 
me — ni < mg < min{n2, ttis) + min(n3, me) 

(11) 

and 



n7 = 715 + 716 — ng 

ris = n2 + 713 - ng 

7714 = 7li + 713 — me 

7717 = 7715 + 7716 — mg 

ms = n2 + 713 - 7719 



(12) 



3.2. FOUR TETRAHEDRA CELL 

As the final ingredient in the K-variable solu- 
tion to the quasi-local updates, we now consider 
four tetrahcdra glued together as shown in fig. 5. 
In the two tetrahedra case we finally had 15 K- 
variables of which 9 were independent and 6 de- 
pendent. The K-variable assignments are sum- 
marised below: 



{OA,AE,OE) 
{DA, AE, DE) 
{OA, AD, OD) 
{OA, AC, OC) 

{DA, AC, DC) 
{BA, AC, BC) 
{BA, AE, BE) 
{OA,AB,OB) 



(ril, n2, n3) 
(7i4, 7i5, n6) 

(7l7, 7l8, 7l9) 

(pl,p2,p3) 

(p4,p5,p6) 
((j4, gS, q&) 
{mA, m5, m6) 
(?7i7, m8, m9) 



(13) 



Of the 24 K-variables 6 are unrestricted, 7 are 



Figure 5. A four-tetrahedra cell. 



restricted but independent and 11 are dependent. 
These are summarised in the following: 



714 



ni + ns - 7i6 



m + n3 > 
ni + 713 > 
"2 + n3 > 



716 > 
7716 > 
P2 > 



(14) 
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min{n2,m5) + ■min{n3,mQ) > mg > uiq — ni 

min{n5 + ne,n2 + ns + pi " P2) > pe > 
min{m5 + me,n2 + ns + pi - P2) > qe > 

(15) 

The new feature at this level is that range re- 
strictions are now of two types: those that de- 
pend only on the unrestricted set as in the case of 
(ng, me,p2) and those that depend on the already 
restricted set as in the case of (rig, mg,pg, gg)- 
The dependent variables are given by: 



714 


= n-i - 


tns - 


ne 


717 


= n5 - 


\-n6 - 


ng 


ns 


= n2 - 


\-n3- 


ng 


rrii 


= Til -1 


- n3 - 


7776 


mj 


= 7715 ^ 


- "76 - 


- mg 


ms 


= 7l2 4 


- "3 - 


mg 


P3 


= 7l2 - 


f 773 - 


P2 


Pi 


= Pl - 


f P3 - 


P6 


P5 


= 775 - 


f - 


P6 


94 


= Pl - 


f P3 - 


96 


95 


= 7775 - 


f me - 


- 96 



4. GLOBEVEN AND GLOBODD 

Finally we are in a position to solve the 
globeven and globodd problems. For globeven, 
one takes two copies of the 4-tetrahedra cell and 
glues them such that the rectangular bases of the 
two cells coincide. As before the number of K- 
variables required will be twice the number for 
each 4-tetrahedra cell minus the number associ- 
ated with the rectangular base. But because of 
the 8 links that are shared between the two cells, 
the independent variables will be reduced by 8. 
The results are similar to eqns(15,16) but more 
lengthy and will not be displayed here. 

For the globodd case we have to take six copies 
of the 4-tetrahedra cells and glue them together 
such that the six rectangular bases now form the 
skin of the 8-cube cluster shown in fig. 2. Now the 
total number of independent K-variables will be 
six times those of the 4-tetrahedra cell reduced by 



the 12 new relations that arise due to the sharing 
of the 12 b-links. 

This completes the solution to the problem of 
finding spin-configurations for the globoeven and 
globodd updates that will respect all the trian- 
gle inequalities. Yet, the resulting equations are 
quite involved and may not result in efficient al- 
gorithms for numerical simulations. A practical 
approach that can be followed to accomplish a Z2- 
flip is to simply increase all the interior spins by 
i except those that are already at the maximum 
value permitted by the code; for these, the spin 
value can be reduced by ^. It is easily seen that 
this "brute force" Z2-Rip will also maintain the 
triangle inequalities. But such moves do not nec- 
essarily respect the balance conditions required 
of Monte Carlo simulations. As a way out of 
this, we swept through the relevant clusters many 
times randomly setting the links to new values al- 
lowed by triangle inequalities. This was achieved 
through the local update methods mentioned be- 
fore. New values were accepted according to a 
metropolis algorithm or a heat bath method car- 
ried over a small chosen subset. The random 
numbers were generated uniformly in the range 
(0, 1).. Whether this method is picking new link 
samples in an unbiased way is yet to be under- 
stood satisfactorily ||]. 
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